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QQ ' 1. Introduction. In this paper, we construct asymptotic expan- 

C"~~- ■ sions for the tail area G of a compound sum, when the summands 

^J . belong to a class of light-tailed subexponential distributions. To be 

(^ I more precise, let X^, i ^ 1, be a sequence of independent random 

^O ■ variables, all having the same distribution F. For any positive in- 



teger n the partial sums Sn = Xi + ■ ■ ■ + X^ have distribution the 
'^ n-fold convolution F*". We set 5o = and therefore F**^ is defined 

_i , as the distribution of the point mass at the origin. Let A'^ be a non- 

H ' negative integer- valued random variable, independent of the X^'s. 

*> • We consider the distribution G of the compound sum Sn, that is 

^ , EF* . Its tail area is G = 1 — G. First order asymptotic results for 

H I G have been obtained by Embrechts, Goldie and Veraverbeke (1979), 

Cline (1987), Embrechts (1985) and Griibel (1987). A second order 
formula may be found in Omey and Willekens (1987). 

Compound sums or subordinated distributions arise as distri- 
bution of interest in several stochastic models. In insurance risk 
theory, it models the total claim amount. For a discussion of is- 
sues related to random sums and insurance risk, we refer to Em- 
brechts, Kliippelberg and Mikosch (1997), Asmussen (1997), Goldie 
and Kliippelberg (1998). Compound sums also appear in queue- 
ing theory, in connection with the stationary distribution of waiting 
times in the GI/G/1 queue. The connection here is not as direct 



as in the insurance risk model in that it is derived from an analysis 
of ladder heights for transient random walks; see, for example, As- 
mussen (1987, p.80). Feller (1971, p.396) and Pakes (1975). Another 
common way in which this model occurs is through the solution of 
a transient renewal equation. An example of this occurs in branch- 
ing processes, where we obtain a geometric-compound sum in the 
analysis of the mean number of particles alive at a given time in 
an age-dependent subcritical process; see Athreya and Ney (1972, 
p. 151). We refer to Feller (1971, chapter XI) for a discussion of 
transient renewal theory. For further applications of subexponen- 
tiality in transient renewal theory, we refer to Teugels (1975) and 
Embrechts and Goldie (1982). 

Throughout the paper, we assume that the XiS are nonnegative. 

2. Main results. If it exists, the hazard rate h = F' / F yields the 
representation of the distribution function F as 



F(t) = F(to)exp('- / h{u)du 



We write Id the identity function on M; for any positive real number 
r, the function Id*^ maps t to V . From the representation of F with 
its hazard rate, we see that if /i ~ a/Id at infinity, then F is regularly 
varying with index —a. If limj^oo h{t) = a then F{t) = e""*^"'^"'""*-^-'-' 
has a tail behavior close to that of an exponential distribution. 
Since we are interested in light subexponential tails, it is natural 
to consider hazard rates such that 

h is regularly varying, 

limt^oo th{t) = -|-oo and Ivait^oo h{t) = . 

In order to be not too close to the Pareto type distributions, we will 
strengthen this assumption by requiring that 

liminft/i(t)/logt > 0. (2.2) 

t — >oc 

This excludes distributions with tail 6"^^°^*-' with a < 2, but 
include those for which a ^ 2. It also includes the subexponential 
Weibull distributions, or more generally, those with tail of the form 
t^e~^ with a positive and less than 1. 



As observed in Barbe and McCorniick (2004, 2005), smoothness 
is a key requirement to obtain asymptotic expansions. For our pur- 
poses, a good class of regularly varying functions are the smoothly 
varying ones of given order, whose definition we now recall. 

Definition. A function h is smoothly varying of index a and order 
m if it is ultimately m-times continuously difjerentiable and its m-th 
derivative is regularly varying of index a — m. 

Clearly, if the hazard rate is m times differentiable, the tail 
function F can be differentiated ttt. + 1 times. 

The next notation we need to introduce pertains to the Laplace 
characters. We write D the derivation operator; that is, if g is 
differentiable, T>g is its derivative. As is customary, we define D*' 
to be the identity, and for any positive integer i we define D* by 
induction as DD*~^. 

We write ^F,i the i-th moment of F. 

Definition. (Barbe and McCormick, 2004)- Let F be a distribution 
function having at least m moments. Its Laplace character of order 
m is the differential operator 

'-F,m — 2_^ ^ ^F.i-LJ • 



Laplace characters have useful algebraic properties which are 
described in Barbe and McCormick (2004). In particular, consider 
the ring ]Rj„[^] defined as the quotient ring of polynomials in 
D modulo the ideal generated by D"*"*"^. Laplace characters are 
elements of this ring, and can be multiplied. It may be helpful to 
think of a Laplace character as a formal Laplace transform Ee~^^ 
where X has distribution F, expressed as a formal Taylor series in D, 
dropping all terms in D'^+i^ j^f^+s^ _ _ _ Then, the multiplication in 
the ring IRm.[L) ] amounts to the usual multiplication of Taylor series, 
dropping any term in 0'"+-'^, 0*"+^, ... In particular, in Mm[E)], we 
have LH*K,m = LH,mLK,m- In what follows, we always consider 
Laplace characters of order m as members of M,ra[D], and all the 
operations on Laplace characters are in that quotient ring. 



The following theorem provides an asymptotic expansion for the 
tail of G. 

Theorem 2.1. Let F be a distribution function whose hazard rate 
is smoothly varying with negative index at least —1 and positive order 
m. Assume further that (2.2) holds and that the moment generating 
function of N is finite in a neighborhood of the origin. Then for any 
nonnegative integer k at most m 

G = ENLp.(iv-i) ,fcF + o(/i^F) . 

Remark. It is shown in Barbe and McCormick (2005, Lemma 
4.1.1) that under the assumptions of Theorem 2.1 the asymptotic 
equivalence F^^^ -^ {—l)^h'^F holds. Therefore, the remainder term 
in the above formula could be written as o{F^^'). 

3. Examples. We illustrate the use of Theorem 2.1, considering 
the cases where N has a Poisson and a geometric distribution. 

Example 1. Assume that N has a Poisson distribution with param- 
eter a. Sums with a Poisson number of summands are commonly 
used in insurance mathematics, modelling total claim size (see Beir- 
lant et al., 1996, Embrechts et al., 1997, Willmot and Lin, 2000). 
The following expansion is easily derived. 

Proposition 3.1. Let F be a distribution function satisfying 

the assumptions of Theorem 2.1. If N has a Poisson distribution 
with parameter a, then G = aLo^rnF + o{h"^F). Moreover, Lo^rn = 

Proof. Combine Theorem 2.1 and the proof of Corollary 4.4.2 in 
Barbe and McCormick (2004) to obtain the expansion aLcmF. To 
obtain the expression for Lcm, write, in the quotient ring, 

a" 






-F,7n — ^ Z^ \^| F,rn. 



The above formula is easily implemented with a computer alge- 
bra system. For example, the following Maple code calculates 

Qga(/.F,„ — Id) 



mu[0] :=1: 

LF:=sum(' (-l)'"j*mu[j] *D'j/j ! ' ,'j'=0. .m+l) : 

taylor(a*exp(a*(LF-l)) ,D=0,m+l) ; 

Setting ?Ti = 3 in the previous code yields the first four terms, 

E{NL^-^) = aid - oVkiD + —{a^i%^ + ^iF,2)^'' 

q2 

- — (O^/^F 1 + 3a;Uj7 i^F,2 + MF,3)D^ • 
D 

To give a very concrete example, assume that F is the Weibull 
distribution with parameter 1/3, so that F{t) = e~* . Define 
er{t) = fe~^ . We obtain, after evaluation of E{NLp^^), and 
using a computer algebra package, 

G = aeo + 2a^e_2/3 + 20^(20 + a)e_4/3 + 4a^(20 + a)e_5/3 

40^(1680 + 60a + a2) ^ o/.^n^ .r. 2^ 

^ -e_2 + Sa^ (1680 + 60a + a^)e_7/s 



3 

2a2 (403200 + 9120a + 140a2 + a^ 



■e-8/3 + o(e_8/3)- 



Perhaps the only remarkable feature of such 7 terms expansion is 
that it can be done. 

Example 2. Motivated by applications to queueing theory (see e.g., 
Cohen, 1972, or Bingham, Goldie and Teugels, 1987, p. 387), consider 
the case where N has a geometric distribution with parameter a, 
that is A'^ is a nonnegative integer n with probability (1 — a)a'^. 
Again, Theorem 2.1 provides a compact expression of the asymptotic 
expansion of G, and the issue is how to actually compute it. 

Any polynomial in D with nonvanishing constant term is invert- 
ible in the quotient ring Mj„[D]. Therefore, 

ENLp^.^-r^^ra = {l-a)Y^ a^nL^p-^ = a(l - a)(Id - aLF,„.)-^ . 

Consequently, the following result holds. 

Proposition 3.2. Let F be a distribution function satistfying the 
assumptions of Theorem 2.1. If N has a geometric distribution with 
parameter a, then G = a(l — a)(Id — aLp^m)~'^F + oiJi^F). 



Setting m = 3, we obtain, as in the previous example, with the 
help of a computer algebra package, with 6 = a/(l — a), 

ENL^-^ = 6Id - 26Vf,iD + fo^(/UF,2 + 36^|. 1)0^ 
52 

- y (126Vf,1 + 9bfiF,ll^F,2 + /UF,3)D3 . 

For instance, when F is the Weibull distribution with parame- 
ter 1/3, the calculation of ENL^^^ yields the following 7 terms 
expansion — expressed solely with a, the formula contains alternat- 
ing signs; it is numerically slightly more stable when expressed with 
h = a/{l-a). 



G = heo + ^b^ e_2/3 + 46^(20 + 36) 6.4/3 + 862(20 + 36) 6.5/3 

^2/1^^ I 1 CJ, I !,2\ „ I inoJ,2/'i^n i ^ r.u i ^,2 



+ 326^(140 + 156 + 6^) e.2 + 1926^(140 + 156 + 6^) 6.7/3 

8/3 + 0(6-8/3) 



+ 806^(6720 + 4566 + 286^ + 6^) 6.8/3 + 0(6-8/3) • 



4. Proof. When m vanishes, Theorem 2.1 is due to Embrechts, 
Goldie and Veraverbeke (1979, p. 342). Therefore, we will prove it 
when m is at least 1. 

It is convenient to introduce a pseudo-semi-norm on tails. If K is 
a distribution function, we write 

li?l^ = sup(ir/F)(t), 

with the convention 0/0 = 0. This generates balls B{F, r) containing 
all tails K which are less than rF. We write B{F) the union of all 
these balls for all positive r. 

We write G„ for the n-fold convolution F*". 

We start by recalling Kesten's global bound on tail function of 
self-convolutions of subexponential distributions; see Athreya and 
Ney (1972, §IV.4, Lemma 7). It asserts that for any positive e there 
exists a positive A such that for all positive integers n, 

|G„|^^^(l + e)". (4.1) 

We also need a precise estimate of the order of magnitude of 
derivatives of F. As noted in the Remark following Theorem 2.1, 
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Lemma 4.1.1 in Barbe and McCormick (2005) shows that for any 
nonnegative k at most m, 

i?W ^(_l)fc/jfci?. (4.2) 

Finally, we also need a basic representation of convolution in terms 
of operators. For any distribution function K with support in the 
nonnegative half- line and any rj positive and less than 1, define the 
operator 



TK,^f{t) = r fit - x) dK{x) . 
Jo 



For any positive c we also define the multiplication operator M^ 
acting on functions by 

MJ{t) = fit/c) . 

These two operators allow us to write a convolution in a way 
suitable for our analysis. Define the powers T^ by 7"^ = Id and 
T^~t = Tk^t] TK,r]- Using Proposition 5.1.1 in Barbe and McCormick 
(2004) inductively we obtain a representation for the distribution 
function, valid on the nonnegative half line, 

(4.3) 
Our first lemma is a simple moment bound. 



Lemma 4.1. Let i he a nonnegative integer, and let e he a positive 
real number. There exists ti such that for any t at least f i and any 
distribution function K in B{F), 



x^dK{x) ^ {l + e)\K\pfF{t)^ 



Proof. For any nonnegative integer z, an integration by parts yields 

/"OO 

X* dK{x) = fK{t) + i / x'-^K{x) dx . (4.4) 



The right hand side of this equality is less than \K\p times the same 
expression with K replaced by F. Consequently, it suffices to prove 
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the result when K is F. In that case, let M be a positive real number 
so that e{M — i) ^ i. Since Id/i tends to infinity at infinity, h is more 
than M/Id ultimately. For any t large enough and any x at least t, 



This implies that the integral in the right hand side of (4.4), when 
F is substituted for K, is at most et^F{t). ■ 

Our next lemma contains the main argument of the proof, namely 
that a Tk,-^ operator is in some sense very close to a Laplace 
character as far as tail behavior is concerned when applied to F 
and its derivatives. 



Lemma 4.2. For any fixed integer p at most m, 

lim sup -'-^= — ■ — = (t) = . 



Proof. The proof of Lemma 4.2.3 in Barbe and McCormick (2005) 
shows that for any 5 positive, 

5/h{t) 

F^P^ {t - x) dK{x) - LK,m-pF^''^ (4.5) 



is at most 

poo 

V |F(P+^')(i)| / x^dK{x) (4.6) 

,-^__^ Js/h(t) 



O^j^m—p 



f'5/h(t) I'X rn-p-1 

+ / / T^ -|F(-)(t-z+y)-F(™)(t)|dydK(x).(4.7) 

Jo Jo \1TT' ~ P ~ ^)- 

Let e be a positive number. Using Lemma 4.1 and (4.2), we see that 
for large t, the term (4.6) is less than 

^('>^w-(/^)'''(4^ 

Since F is rapidly varying, this is ultimately less than e\K\ph'^F. 
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The proof of Lemma 4.2.3 in Barbe and McCormick (2005) shows 
that for 5 small enough, for any t large enough and for any K 
in B{F), the double integral (4.7) is at most €\K\piiF,m-ph"^F. 
Hence, we have shown that (4.5) is at most e\K\p{nF^rn-p + l)h'^F 
ultimately uniformly over B{F). 

The proof of Lemma 4.2.4 in Barbe and McCormick (2005) shows 
that for any positive S and rj, ultimately uniformly over B{F), 



l5/h{t) 

This proves Lemma 4.2 



f \F^P\t-x)\dK{x)^e\K\ph"'F{t). 
Js/hd) 



Lemma 4.2 yields the following estimate on an operator T com- 
posed with a Laplace character applied to a derivative of F. 

Lemma 4.3. The following uniform limit holds: 

1. ' K,rii-H,m — p^ '-K*H,m — pJ- „ 

lim sup = = = . 

t^oo , \K\ h"^F V ^"-^ 



Proof. Since TV, 77 is linear and 

O^j^m—p 

the result follows from Lemma 4.2 and Lemma 2.1.4 in Barbe and 
McCormick (2004). ■ 

The next two lemmas will take care of some remainder terms. 
The first one asserts that terms of order o{h"^F) remain so through 
the action of some T operators. 



Lemma 4.4. Let q he a nonnegative integer and e be a positive 
real number. There exist t2, some positive A and rj, such that for 
any positive integer i, 

rp^^{h''F)^A{l + eyh'iF 
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on [^2, oo)- 

Proof. Let e be a positive real number. Since h is regularly varying 
with negative index, provided ry is small enough, h{t—x) ^ (l+e)/i(t) 
for any t large enough and any x nonnegative and at most rjt. 
Therefore, for t at least t'2, 



rv^ 

TF,r,{h'^F){t) = / /i«F(t - x) dF(x) 

^ (l + e)/i''(t) I F{t-x)dF{x) 
Jo 

^ il + e)h^t)F*^{t). 

By induction, it follows that 

Ti,^^{h^F){t) ^ (l + e)^(/i'^F*(^+i))(t). 

Using Kesten's bound, (4.1) above, this yields that Tp.^{h'^F) is 
ultimately at most A{1 + e)^*/i'^F, finishing the proof since e is 
arbitrary. ■ 

Our penultimate lemma will be used to handle the terms involving 
the multiplication operators in (4.3). 

Lemma 4.5. Let e be a positive real number. There exists t^ such 
that for any positive integers i and m, 

|/Wi/^FMi/(i_^)G,| ^ (1 + eyh^'+'F 

on [^3, 00). 

Proof. Kesten's bound shows that 

\F{t7])G,{t{l -v))\^ F{t7])A{l + eyF{t{l - r/)) . 

Arguing as in Lemma 4.2.1 in Barbe and McCormick (2005), 
F{tri)F(t{l ~ T])) is o{h^F{tyj for any positive q. This implies the 
result. ■ 

Our last lemma is stated merely to avoid digression in the argu- 
ment later on. 
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Lemma 4.6. Let e be a positive number. There exists A such that 
for any positive integer n 

Proof. The lemma folllows from Marcinkiewicz-Zygmund's inequal- 
ity (see Chow and Teicher, 1988, §10.3, Theorem 3), which implies 
that /iG„,j ^ An^ for some constant A. ■ 

We can now conclude the proof of Theorem 2.1. Combining 
Lemmas 4.4 and 4.5, there exists an interval [^3,00) on which for 
any j and k with ^ j ^ k ^ n, any positive i and n with i ^ n, 

I^f:^'(Wi/^^Wi/(i-^)G„-0| ^ ^(1 + e)"/i'"+^F. 

Representations (4.3) yield, on [t3,oo), 

\Gn- E njTG„_^,,_T\^Anil + erh"^+'F. (4.8) 

Let e be a positive real number, small enough so that E{1 + 2e)^ 
is finite. Let 5 be a positive real number. Combining Lemmas 4.2, 
4.3, 4.4 and 4.6, using also Kesten's bound, ultimately, uniformly in 
n and i at most n, 

I 7-i— 1 7- 7? T^~^ I 7?l 

I ' F,J7 'G„-i,l-ri^ ' F,ri '-G„-i+i,m^ \ 

^ "'X'? \i'^Gr.^,,l-n-LG„_,,m)F\+Vf^^ |(7F,,,/-G„_,,m-/-G„_,+i,„)^| 

i^25{A^ + A){l + e)''h'^F . (4.9) 

Using the same combination of lemmas, we also have, ultimately, 
uniformly in n and j at most n — 1, 

\n,,LG„-,-umF - n:,'l-G^-,,^F\ ^ A\l + er5h^F. (4.10) 

We take A to be at least 1, simply to ensure that A^ is more than 
A. Summing (4.10) for j positive and less than i and adding (4.9), 
we obtain 

I n:^ Tg„_^,i-vF - LG^_„mF\ ^ AAHh^Ft{l + e)" 
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on some interval [t4,oo). Summing these inequalities for i positive 
and at most n and combining with (4.8) yield 

\Gn-nLG.^_^^^F\^WA^6n{n + l){l + erh^F. 

Since the moment generating function of N is finite at log(l + 2e) 
and S is arbitrary, Theorem 2.1 follows. ■ 
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